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QO ■ Abstract 

We study a particular class of relativistic nuclear energy density functionals in which only nucleon 
degrees of freedom are explicitly used in the construction of effective interaction terms. Short- 
^ ' distance (high-momentum) correlations, as well as intermediate and long-range dynamics, are 

encoded in the medium (nucleon density) dependence of the strength functionals of an effective 
, interaction Lagrangian. Guided by the density dependence of microscopic nucleon self-energies 

I in nuclear matter, a phenomenological ansatz for the density-dependent coupling functionals is 

m ' 



X 



accurately determined in self-consistent mean-field calculations of binding energies of a large set 
of axially deformed nuclei. The relationship between the nuclear matter volume, surface and 



o 

00 . 

, symmetry energies, and the corresponding predictions for nuclear masses is analyzed in detail 



The resulting best-fit parametrization of the nuclear energy density functional is further tested 
in calculations of properties of spherical and deformed medium-heavy and heavy nuclei, including 
binding energies, charge radii, deformation parameters, neutron skin thickness, and excitation 
energies of giant multipole resonances. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.10.Dr, 21.10.Ft 
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I. INTRODUCTION 



Among the microscopic approaches to the nuclear many-body problem, the framework 
of nuclear energy density functionals (NEDF) provides the most complete and accurate 
description of ground-state properties and collective excitations over the whole nuclide chart. 
Probably no other method achieves comparable accuracy at the same computational cost. 
At the level of practical applications the NEDF framework is realized in terms of self- 
consistent mean-field (SCMF) models based, for instance, on the Gogny effective interaction, 
the Skyrme energy functional, and the relativistic meson-exchange effective Lagrangian 
[21. In the mean-field approximation the dynamics of the nuclear many-body system is 
represented by independent nucleons moving in self-consistent potentials, that correspond 
to the actual density and current distributions of a given nucleus. 

The SCMF approach to nuclear structure is analogous to Kohn-Sham density functional 
theory 3|, y], and nuclear mean- field models approximate the exact energy functional, which 
includes all higher-order correlations, with powers and gradients of ground-state nucleon 
densities and currents j^. In particular, a number of very successful relativistic mean-field 
(RMF) models have been constructed based on the framework of quantum hadrodynamics 
(QHD) . There are important advantages in using functionals with manifest covariance 
8j. The most obvious is the natural inclusion of the nucleon spin degree of freedom, and 
the resulting nuclear spin-orbit potential which emerges automatically with the empirical 
strength in a covariant formulation. The consistent treatment of large, isoscalar, Lorentz 
scalar and vector self-energies provides a unique parametrization of time-odd components 
of the nuclear mean-field, i.e. nucleon currents, which is absent in the non-relativistic rep- 
resentation of the energy density functional. The empirical pseudospin symmetry in nuclear 
spectroscopy finds a natural explanation in terms of relativistic mean fields jsj. On a micro- 
scopic level, it has been argued jsl that a covariant formulation of nuclear dynamics manifests 
the true energy scales of QCD in nuclei, and is consistent with the nonlinear realization of 
chiral symmetry through the implicit inclusion of pion-nucleon dynamics in the effective nu- 
cleon self-energies. A covariant treatment of nuclear matter provides a distinction between 
scalar and four-vector nucleon self energies, leading to a very natural saturation mechanism. 

In conventional QHD the nucleus is described as a system of Dirac nucleons coupled 
to exchange mesons through an effective Lagrangian. The isoscalar scalar a meson, the 



2 



isoscalar vector uj meson, and the isovector vector p meson build the minimal set of meson 
fields that is necessary for a description of bulk and single-particle nuclear properties. Of 
course, at the scale of low-energy nuclear structure, heavy-meson exchange is just a con- 
venient representation of the effective nuclear interaction. At the energy and momentum 
scales characteristic of nuclei, the only degrees of freedom that have to be taken into account 
explicitly in the description of many-body dynamics are pions and nucleons. The behavior 
of the nucleon-nucleon (NN) interaction at long and intermediate distances is determined 
by one- and two-pion exchange processes. The exchange of heavy mesons is associated with 
short- distance dynamics that cannot be resolved at low energies that characterize nuclear 
binding and, therefore, can be represented by local four-point (contact) NN interactions, 
with low-energy (medium-dependent) parameters adjusted to nuclear data. These con- 
cepts of effective field theory and density functional theory methods have recently been 
used to derive a microscopic relativistic energy density functional framework constrained 
by in-medium QCD sum rules and chiral symmetry 10, lUj. The density dependence of 
the effective nucleon-nucleon couplings is determined from the long- and intermediate-range 
interactions generated by one- and two-pion exchange processes. They are computed us- 
ing in-medium chiral perturbation theory, explicitly including A(1232) degrees of freedom 



12] . Regularization dependent contributions to the energy density of nuclear matter, calcu- 



lated at three-loop level, are absorbed in contact interactions with parameters representing 
unresolved short-distance dynamics. 

However, even in a fully microscopic approach that starts from a description of symmetric 
and asymmetric, homogeneous and inhomogeneous nuclear matter, the parameters of a 
nuclear energy density functional have still to be fine tuned to structure data of finite 
nuclei. This is simply because gross properties of infinite nuclear matter cannot determine 
the density functional on the level of accuracy that is needed for a quantitative description 
of structure phenomena in finite nuclei. For most functionals this tuning is performed 
on a relatively small set of spherical closed-shell nuclei, mainly because they are simple 
to calculate and can therefore be easily included in multiparameter least-squares fits. A 
problem arises, however, because ground-state data of closed-shell nuclei include long-range 
correlations that cannot really be absorbed into mean-field functionals. Generally this will 
affect the predictive power of energy density functionals when they are used in the description 
of phenomena related to the evolution of shell structure. For instance, soft potential energy 
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surfaces and/or small energy differences between coexisting minima in deformed nuclei, 
are often difficult to describe using functionals adjusted solely to data of spherical nuclei, 
even when sophisticated models are employed that include angular momentum and particle 
number projection, as well as intrinsic configuration mixing. 

In this work we explore a class of relativistic energy density functionals originally intro- 



duced in Refs. 
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IfJ but, instead of using low-energy QCD constraints for the medium 



dependence of the parameters, a phenomenological ansatz is adjusted exclusively to masses 
of a relatively large set of axially deformed nuclei. The phenomenological approach, although 
guided by microscopic nucleon self-energies in nuclear matter, gives us more freedom to in- 
vestigate in detail the relationship between global properties of a nuclear matter equation of 
state (volume, surface, and asymmetry energies) and the corresponding predictions for nu- 
clear binding energies. Eventually the goal will be to develop an energy density functional 
that does not implicitly contain symmetry breaking corrections and quadrupole fluctua- 
tion correlations, and is therefore better suited for the new relativistic model that uses the 
generator coordinate method to perform configuration mixing of angular-momentum and 

nn 

particle-number projected relativistic wave functions pJl, IMj. The idea is that those corre- 
lations that we wish to treat explicitly, should not be included in the density functional in an 
implicit way, i.e. by adjusting parameters to data which already include correlations. The 
solution could be to adjust the functional to pseudodata, obtained by subtracting correla- 
tion effects from experimental masses and, eventually, radii. This is most easily done using 
masses of axially deformed nuclei with large deformation parameters, because the dominant 
contribution to their ground-state correlation energies is the rotational energy correction 
[isl . which is relatively simple to calculate. Approximate methods for the calculation of 
correlations have been developed [igI , that enable a systematic evaluation of correlation en- 
ergies for the nuclear mass table. Of course one expects that the corresponding modifications 
of the parameters of the energy density functional will be relatively small, but even a small 
change in the relative contribution of various interaction terms could be the decisive factor 
in the description of soft potential energy surfaces, coexistence of prolate and oblate shapes, 
level ordering, etc. Very recent examples include the phenomenon of shape coexistence in 



neutron-deficient Kr isotopes [l7|, [l8| , and the description of singular properties of excitation 



spectra and transition rates at critical points of quantum shape phase transitions [19]. As 
a first step toward the construction of a relativistic density functional that could provide a 
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more accurate description of phenomena related to the evolution of shell structure, in this 
work we explore the possibility to determine the parameters of a given functional using only 
binding energies of axially deformed nuclei. 

In Section [TTl we construct the relativistic nuclear energy density functional based on 
the above conjectures, and discuss the necessary approximations and fitting strategies. In 
Section llllt starting from microscopic nucleon self-energies in nuclear matter, and empir- 
ical global properties of the nuclear matter equation of state, the functional is accurately 
determined in a careful comparison of the predicted binding energies with data, for a set 
of 64 axially deformed nuclei in the mass regions A ^ 150 — 180 and A ^ 230 — 250. In 
Section IIVI the new energy density functional is thoroughly tested in a series of illustrative 
calculations of properties of spherical and deformed medium-heavy and heavy nuclei, in- 
cluding binding energies, charge radii, deformation parameters, neutron skin thickness, and 
excitation energies of giant multipole resonances. Section |V] summarizes the results of the 
present investigations and ends with an outlook for future studies. 



II. RELATIVISTIC NUCLEAR ENERGY DENSITY FUNCTIONAL 



The basic building blocks of a relativistic nuclear energy density functional are the den- 
sities and currents bilinear in the Dirac spinor field ip of the nucleon: 

i>OrTi> , Ore{l,Ti}, r G {1,7^,75,757^,(7^^} . (1) 

Here Xj are the isospin Pauli matrices and F generically denotes the Dirac matrices. The 
nuclear ground-state density and energy are determined by the self-consistent solution of 
relativistic linear single-nucleon Kohn-Sham equations. To derive those equations it is useful 
to construct an interaction Lagrangian with four-fermion (contact) interaction terms in the 
various isospace-space channels: 

isoscalar-scalar: {ipipY 

isoscalar- vector: {ilj'y^ip){ip'y'^ip) 

isovector-scalar: {ipTip) ■ {ipTip) 

isovector-vector: [ipf'^^^il)) ■ {ipT'-f'^ilj) . 
Vectors in isospin space are denoted by arrows. A general Lagrangian can be written as a 
power series in the currents ipOrTip and their derivatives, with higher-order terms repre- 



senting in-medium many-body correlations [20 
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24| . We will adopt the approach 



of Refs. [lO|, lull and construct a Lagrangian with second-order interaction terms only, with 
many-body correlations encoded in density- dependent coupling functions. In complete anal- 
ogy to the successful meson-exchange RMF phenomenology, in which the isoscalar scalar a 
meson, the isoscalar vector lo meson, and the isovector vector p meson build the minimal set 
of meson fields that is necessary for a description of bulk and single-particle nuclear proper- 
ties, we consider an effective Lagrangian that includes the isoscalar-scalar, isoscalar-vector 
and isovector-vector four-fermion interactions: 

L = ■ d — m)ip 

- ^5s{d,i^i^Wij^) - eV^7 ■ ^^^^^ • (2) 

In addition to the free-nucleon Lagrangian and the point-coupling interaction terms, when 
applied to nuclei, the model must include the coupling of the protons to the electromagnetic 
field. The derivative term in Eq. ([2]) accounts for leading effects of finite-range interactions 
that are crucial for a quantitative description of nuclear density distribution, e.g. nuclear 
radii. Similar interactions can be included in each space-isospace channel, but in practice 
data on charge radii constrain only a single derivative term, for instance 5s{di,ipip){d^ipil)). 
The coupling parameter 5s has been estimated, for instance, from in-medium chiral pertur- 



bation calculation of inhomogeneous nuclear matter [12|]. In the region of nucleon densities 
relevant for the description of finite nuclei (0.1 fm~^ < p < 0.15 fm^^), the coupling strength 
of the derivative term displays a rather weak density dependence and can be approximated 
by a constant value 5s between —0.85 fm^ and —0.7 fm^. Note that the inclusion of an ad- 
justable derivative term only in the isoscalar-scalar channel is consistent with conventional 
meson-exchange RMF models, in which the mass of the fictitious a meson is adjusted to 
nuclear matter and ground-state properties of finite nuclei, whereas free values are used for 
the masses of the uj and p mesons. 

The point-coupling Lagrangian Eq. ([2]) does not include isovector-scalar terms. In the 
meson-exchange picture this channel is represented by the exchange of an effective 5 me- 
son, and its inclusion introduces a proton-neutron effective mass splitting and enhances the 
isovector spin-orbit potential. However, in calculations of ground-state properties of finite 



nuclei, using both meson-exchange 25|, |26| and point-coupling 2J] models, it has not been 



possible to constrain the parameters of the effective interaction in the isovector-scalar chan- 



nel. Although the isovector strength has a relatively well-defined value, the distribution 
between the scalar and vector channels is not determined by ground-state data. To reduce 
the number of adjustable parameters, the isovector-scalar channel may be omitted from 
an energy density functional that will primarily be used for the description of low-energy 
nuclear structure. 

The strength parameters of the interaction terms in Eq. ([2]) are functions of the nucleon 
4-current: 

f = ^^^^ = pu^ , (3) 

where is the 4-velocity defined as (1 — v'^)^^^'^(l,v). In the rest-frame of the nuclear 
system: v = 0. The single- nucleon Dirac equation, the relativistic analogue of the Kohn- 
Sham equation, is obtained from the variation of the Lagrangian with respect to ip: 

Mzd^ - - S^) - (m + Es)] = , (4) 

with the nucleon self-energies defined by the following relations: 

= avip.)f + e^^^^A^ (5) 

^«-2pa dp^'^ dp^"^ + dp / 

= ois{Pv)Ps - ^^Dps (7) 
S^^ = aTv{Pv)f ■ (8) 

In addition to the contributions of the isoscalar-vector four-fermion interaction and the 
electromagnetic interaction, the isoscalar-vector self-energy S'^ includes the "rearrangement" 
terms S^, arising from the variation of the vertex functionals as-, ay, and arv with respect 
to the nucleon fields in the density operator p. The inclusion of the rearrangement self-energy 
is essential for energy-momentum conservation and the thermodynamical consistency of the 
model 26|, |27|, |28(] . and denote the isoscalar-scalar and isovector- vector self-energies, 
respectively. 

In the relativistic density functional framework the nuclear ground state |(^o) is repre- 
sented by the mean- field self-consistent solution of the system of equations (HI) - (IHl), with 
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the isoscalar and isovector 4-currents and scalar density: 

N 

jfi = (</'o|V'7/.V'l0o) = '^klt^i^k , (9) 



k=l 
N 



i = {<f>o\^7t,r^\(f)o) = X^^fc ^kli^T^pk , (10) 

k=l 

N 

Ps = (0o|#|0o) = ^vltpk^k, (11) 

k=l 

where ipk are Dirac spinors, and the sum runs over occupied positive-energy single-nucleon 
orbitals, including the corresponding occupation factors The single-nucleon Dirac equa- 
tions are solved self-consistently in the "no-sea" approximation that omits the explicit contri- 
bution of negative-energy solutions of the relativistic equations to the densities and currents. 
Vacuum polarization effects are imphcitly included in the adjustable density-dependent pa- 
rameters of the theory. 

A large part of this work will be devoted to adjusting the free parameters of the medium- 
dependent point-coupling functionals as, ay, and utv, and the strength 6s of the derivative 
term. To establish the density dependence of the couplings one could start from a microscopic 
(relativistic) EoS of symmetric and asymmetric nuclear matter, and map the corresponding 
nucleon self-energies on the mean-field self-energies Eqs. ([H]) - ([H]) that determine the single- 
nucleon Dirac equation (jl]). This approach has been adopted, for instance, in RMF models 



based on Dirac-Brueckner-Hartree-Fock self-energies in nuclear matter 
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28 



29|, or on 



in-medium chiral perturbation theory (ChPT) calculations of the nuclear matter EoS 10|, 



Jjj. In general, however, energy density functionals determined directly from a microscopic 
EoS do not provide a very accurate description of data in finite nuclei. The reason, of 
course, is that a calculation of the nuclear matter EoS involves approximation schemes and 
includes adjustable parameters that are not really constrained by nuclear structure data. 
The resulting bulk properties of infinite nuclear matter (saturation density, binding energy, 
compression modulus, asymmetry energy) do not determine uniquely the parameters of 
nuclear energy density functionals, which usually must be further fine-tuned to ground-state 
data (masses and/or charge radii) of spherical nuclei. 

In a phenomenological construction of a relativistic energy density functional one starts 
from an assumed ansatz for the medium dependence of the mean-field nucleon self-energies, 
and adjusts the parameters directly adjusted to data of spherical nuclei. This procedure 



was used, for instance, in the construction of the relativistic density-dependent interactions 
TW-99 I26], DD-MEl fl, DD-ME2 30], PKDD 31 1, PKOl [s^. 

This work adopts a different strategy and determines the parameters of the point-couphng 
Lagrangian Eq. ([2]) exclusively from a large data set of binding energies of deformed nuclei. 
First one notes that calculated masses of finite nuclei are primarily sensitive to the three 
leading terms in the empirical mass formula: volume, surface and symmetry energy 

B.E. = a,A + asA^/^ + a j"^ ~ +■■■ . (12) 

Therefore one can generate families of effective interactions that are characterized by differ- 
ent values of a^,, and 04, and determine which parametrization minimizes the deviation 
from the empirical binding energies of a large set of deformed nuclei. This approach differs 
considerably from the standard procedure of fitting parameters of nonrelativistic Skyrme or 
RMF functionals, in which a given set of parameters is adjusted simultaneously to a favorite 
nuclear matter EoS and to ground-state properties of about 10 spherical closed-shell nuclei. 
Deformed systems have generally not been included in fits of parameters of self-consistent 
RMF models. The reason, of course, is that calculation of deformed nuclei is computa- 
tionally more demanding and requires advanced computer codes. In this work parameters 
of relativistic energy density functionals are for the first time directly adjusted to binding 
energies of axially deformed nuclei in the mass regions Ak, 150 — 180 and A ^ 230 — 250. 

To determine the functional form of the density dependence of the couplings as, ay, and 
aTv, one can start from microscopic nucleon self-energies in nuclear matter. In a recent 



analysis of relativistic nuclear dynamics 



33j, modern high-precision nucleon-nucleon (NN) 



potentials (Argonne Vis, Bonn A, CD-Bonn, Idaho, Nijmegen, Viowk) were mapped on a 
relativistic operator basis, and the corresponding relativistic nucleon self-energies in nuclear 
matter were calculated in Hartree-Fock approximation at tree level. A very interesting re- 
sult is that, at moderate nucleon densities relevant for nuclear structure calculations, all 
potentials yield very similar scalar and vector mean fields of several hundred MeV magni- 
tude, in remarkable agreement with standard RMF phenomenology: at saturation density 
a large and attractive scalar field ^ —400 MeV, and a repulsive vector field S^, ^ 350 
MeV. The different treatment of short-distance dynamics in the various NN potentials leads 
to slightly more pronounced differences between the corresponding self-energies at higher 
nucleon densities. Generally, however, all potentials predict a very similar density depen- 
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dence of the scalar and vector self-energies. In the chiral effective field theory framework, in 
particular, these self-energies are predominantly generated by contact terms that occur at 
next-to-leading order in the chiral expansion. 

Of course at the Hartree-Fock tree level these NN potentials do not yield saturation of 
nuclear matter. Nevertheless, the corresponding self-energies can be used as the starting 
point in the modeling of medium dependence of a relativistic nuclear energy density func- 
tional. Guided by the microscopic density dependence of the vector and scalar self-energies, 
we choose the following practical ansatz for the functional form of the couplings 

a,{p) = a, + (6, + c,a;)e-'^'^ {i = S, V, TV) , (13) 

with X = p/psat, and psat denotes the nucleon density at saturation in symmetric nuclear 
matter. Note that the corresponding self-energies are defined in Eqs. (E]) - (IHl). In the next 
section we will adjust the parameters of the ansatz Eq. f|T3l) simultaneously to infinite and 
semi-infinite nuclear matter, and to binding energies of deformed nuclei. The resulting self- 
energies in nuclear matter will eventually be compared to our starting approximation: the 



Hartree-Fock scalar and vector self-energies of the Idaho N^LO potential [SJ]. 

In the isovector channel the corresponding Hartree-Fock tree level nucleon self-energies, 
obtained by directly mapping microscopic NN potentials on a relativistic operator basis, 
are presently not available. Therefore, as it was done in the case of the finite-range meson 



exchange interactions TW-99 



26|, DD-MEl |27|, DD-ME2 |30|, and PKOl [3^, the density 



dependence of the isovector- vector coupling function is determined from the results of Dirac- 
Brueckner calculations of asymmetric nuclear matter [29]. Accordingly, in Eq. f[T^ for the 
isovector channel we set two parameters to zero: ayy = and ctv = 0, and adjust brv and 
(Itv to empirical properties of asymmetric matter and to nuclear masses, together with the 
parameters of the isoscalar channel. 

III. THE EFFECTIVE DENSITY-DEPENDENT INTERACTION DD-PCl 
A. Infinite and semi-infinite nuclear matter 

The usual procedure in the construction of an effective mean-field interaction is the least- 
squares adjustment of parameters to both nuclear matter EoS and to ground-state data 
(masses, charge radii) of spherical nuclei. Instead we generate sets of effective interactions 
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with different values of tlie volume a^, surface a^, and symmetry energy 04 in nuclear matter, 
and analyze the corresponding binding energies of deformed nuclei with A ^ 150 — 180 
and A ^ 230 — 250. The nuclear matter saturation density, compression modulus, and 
Dirac mass will be kept fixed throughout this analysis. The calculated binding energies of 
finite nuclei are not very sensitive to the nuclear matter saturation density, and we take 
Psat = 0.152 fm~^, in accordance with values predicted by most modern relativistic mean- 
field models. In particular, this value has also been used for the meson-exchange effective 
interactions DD-MEl , and DD-ME2 [3^ . From these interactions we also take the Dirac 
effective nucleon mass m|) = m + S5 = 0.58m. In RMF theory the Dirac mass is closely 
related to the effective spin-orbit single-nucleon potential, and empirical energy spacings 



between spin-orbit partner states in finite nuclei determine a re^ 



allowed values: 0.57 < m*j~,/m < 0.61. In a recent study 35j of the relation between 



atively narrow interval of 



finite-range (meson-exchange) and zero-range (point-coupling) representations of effective 
RMF interactions we have shown that, to reproduce experimental excitation energies of 
isoscalar giant monopole resonances, point-coupling interactions require a nuclear matter 
compression modulus Knm ~ 230 MeV, considerably lower than values typically used for 
finite-range meson-exchange relativistic interactions. Thus we take Knm = 230 MeV for all 
effective interactions considered in the present analysis. 

Of course if only nuclear matter properties at the point of saturation density were speci- 
fied, one could parametrize a number of realistic effective interactions that would be difficult 
to compare at the level of finite nuclei. In particular, nuclear structure data do not con- 
strain the nuclear matter EoS at high nucleon densities. Therefore, in addition to psat, 
m^, and Knm) "we fix two additional points on the E{p) curve in symmetric matter to the 
microscopic EoS of Akmal, Pandharipande and Ravenhall fl, based on the Argonne V,s 
NN potential and the UIX three-nucleon interaction. This EoS has extensively been used 
in studies of high-density nucleon matter and neutron stars. At almost four times nuclear 
matter saturation density, we choose the point p = 0.56 fm~^ with E/A = 34.39 MeV and, 
to have an overall consistency, one point at low density: p = 0.04 fm~^ with E/A = —6.48 
MeV (cf. Table VI of Ref. 36|). As we have already emphasized in the previous section, 
by adjusting mean-field interactions exclusively to a microscopic EoS like, for instance, the 



one calculated in Ref. 



36| , it is not possible to obtain a very accurate description of nuclear 



structure. Ground-state nuclear data must be used to fine tune the parameters of effective 
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interactions. 

In contrast to the Dirac mass and saturation density, the nuclear matter volume energy 
coefficient has a decisive influence on the calculated binding energies of finite nuclei. 
Using the framework of nonrelativistic Skyrme functionals, it was recently shown that even 
a relatively small change in the volume energy 0.5%) can have a pronounced effect on the 
calculated masses of heavy and superheavy nuclei, as compared with experimental values 



15|, 



37l |. In the framework of RMF models no attempt has been made so far to constrain 



the value of volume energy better than the interval —16.2 MeV < a^, < — 16 MeV. To study 
in more detail the effect of volume energy on masses, we have generated point- coupling 
effective interactions characterized by the following values of the coefficient: = —16.02 
MeV (set A), = -16.04 MeV (set B), a„ = -16.06 MeV (set C), a„ = -16.08 MeV 
(set D) , = -16.10 MeV (set E), = -16.12 MeV (set F), a„ = -16.14 MeV (set 
G) and a„ = —16.16 MeV (set H). The corresponding parameters of the ansatz Eq. ( IT3|) 
for the functional form of the isoscalar couplings, are collected in Table [H Note that to 
reduce the number of free parameters, we have set the value cy = 0. The resulting binding 
energy curves for symmetric nuclear matter are plotted in Fig. [1], together with the EoS of 
the meson-exchange effective interaction DD-ME2 , and the microscopic EoS of Ref. 36|. 
The two points on the microscopic EoS that have been used to adjust the parameters are 
represented by large ffiled circle symbols. Because of the anchor at p = 0.56 fm~^, the 
new binding energy curves are, of course, different from DD-ME2 and much closer to the 
microscopic EoS. However, the high-density behavior does not influence much the description 
of low-energy nuclear structure data. 

The isovector channel of the energy density functional determines the density dependence 
of the nuclear matter symmetry energy 



S2ip) = 04 + ^(p - psat) + ^^(P - Psatf + " " " • (W) 



The parameter po characterizes the linear density dependence of the symmetry energy, and 

AKq is the isovector correction to the compression modulus. Experimental masses, unfor- 
;unately, do not place very strict constraints on the parameters of the expansion of 5*2 (p) 
381], but self-consistent mean-field calculations show that binding energies can restrict the 

values of 5*2 at nucleon densities somewhat below saturation density, i.e. at p ~ 0.1 fm~^. 

Additional information on the symmetry energy can be obtained from data on neutron skin 
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thickness and excitation energies of giant dipole resonances. Although values of neutron 
radii are available only for a small number of nuclei and the corresponding uncertainties are 
large, recent studies have shown that relativistic effective interactions with volume asym- 
metry 04 in the range 31 MeV < 04 < 35 MeV predict values for neutron skin thickness 
that are consistent with data, and reproduce experimental excitation energies of isovector 
giant dipole resonances (cf. Ref. 39| and references therein cited). Therefore we keep the 
volume asymmetry fixed at 04 = 33 MeV, and vary the symmetry energy at a density that 
corresponds to an average nucleon density in finite nuclei: (p) = 0.12 fm~^. The quantity 
S2{p = 0.12 fm"^) will be denoted (^2)- 

Calculated binding energies and charge radii are strongly influenced by the choice of the 
surface energy coefficient a^. In the present model the value of this quantity is determined 
by the strength 6s of the derivative coupling term in the point-coupling Lagrangian Eq. ([2]). 
For each effective interaction (sets A-H), we have calculated the surface energy and surface 
thickness of semi- infinite nuclear matter [4^, for several values of the parameter 6s in the 
range —0.76 fm'* < 5s < —0.86 fm'*. In Fig. [2] we plot the corresponding surface energies 
as functions of the surface thickness, in comparison to the point obtained with the finite- 
range interaction DD-ME2 (t = 2.108 fm, = 17.72 MeV). Considering that DD-ME2 
has an rms error of only 0.017 fm when compared to data on absolute charge radii and 



charge isotope shifts 



30], and also taking into account the comparison between DD-ME2 



and point-coupling RMF interactions of Ref. [33] , the following range for the parameter of 



the derivative coupling term can be deduced: -0.80 fm^ > 6s > -0.84 fm^ m very good 
agreement with the microscopic estimate of Ref. 121] for the region of nucleon densities 
p ^ 0.1 fm^. 



B. Deformed nuclei 



If an effective interaction is adjusted to masses of finite nuclei by varying the volume, 
symmetry, and surface energies, the parameters of the energy density functional that deter- 
mine these quantities will generally be correlated because of Eq. f|T2l) . When only a small 
number of nuclei is considered, satisfactory results can be obtained with various, in general 
linearly dependent combinations of parameters. The new effective point-coupling interac- 
tions will therefore be analyzed on a set of 64 deformed nuclei, listed in Table [iTl To resolve 
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the surface and volume contributions to binding energy, nuclides with mass number ranging 
from 154 to 250 are considered. The variation of the asymmetry coefficient 

in the range from 0.018 to 0.054, should suffice to deduce the isovector parameters that 
govern the symmetry energy contribution. The effect of shell closure is minimized by taking 
into account only well deformed nuclei. Pairing correlations are treated in the BCS approxi- 
mation with empirical pairing gaps (5-point formula). The pairing model space includes two 
major oscillator shells (2/za;o) above the Fermi surface. The self-consistent single-nucleon 
RMF equations are solved by expanding nucleon spinors in terms of eigenf unctions of a 



deformed, axially symmetric harmonic oscillator potential [4l|. In this work calculations 
of nuclear ground states are performed in a large basis of 16 major oscillator shells, and 
convergence has been tested in several calculations with 18 oscillator shells. After the self- 
consistent equations are solved, the microscopic estimate for the center-of-mass correction 
is subtracted from the total binding energy 

^- = -^- 

where Pcm is the total momentum of a nucleus with A nucleons. 

For each effective interaction with given volume energy a„ (sets A-H), and for six values 
of the symmetry energy (^2) = 27.6, 27.8, 28.0, 28.2, 28.4, and 28.6 MeV, we have adjusted 
the surface energy (i.e. the coupling strength 63 of the derivative term in the Lagrangian 
Eq. ([2])) to a value that minimizes the deviation of the calculated binding energies from data, 
for the set of nuclei listed in Table HTl The required accuracy is 0.05%, which approximately 
corresponds to an absolute error of ±1 MeV for the total binding energy. The resulting 
surface energies are plotted in Fig. [3] as functions of the volume energy, for each value of 
the symmetry energy (S'2). At this point we have a set of 48 parametrizations of the energy 
density functional. Fig. H] displays the corresponding x^- values 



where E^^(i) denote experimental binding energies 43|], E^^{i) are the corresponding the- 
oretical values, AE^^{i) = 0.0005E'^^{i), and the sum runs over the set of 64 deformed 
nuclei. Although the span of x^-values is very large, the functional dependence of on 
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is smooth and, for each value (5*2 ) of the symmetry energy, there is a unique combination of 
volume and surface energies that minimizes x^- The minima of each curve are collected in 
Fig. O Also in this plot vs displays a smooth parabola, with the absolute minimum at 
the point = -16.06 MeV, (^2) = 27.8 MeV, and a, = 17.498 MeV. The x^-values of the 
neighboring points are not much larger, but obviously the systematics of binding energies 
excludes effective interactions with < —16.10 MeV. 

This result is illustrated in much more detail in Figs. [6lfTT| where we display the absolute 
deviations of the calculated binding energies from the experimental values, for the effective 
interactions that correspond to each of the points included in Fig. [5l Because these inter- 
actions have already been optimized with respect to ag (cf. Fig. [3]) and {S2) (cf. Fig. Sj), 
Figs. IMTT] show the isospin asymmetry (a^), and mass dependences of the absolute errors 
of calculated binding energies as functions of volume energy at saturation a^. Positive devi- 
ations correspond to under-bound nuclei. We notice that the interaction with = —16.06 
MeV (cf. Fig. [7]) not only corresponds to the lowest value, but also that it is the only 
one which does not display any visible isotopic or mass dependence of the deviations of 
calculated masses. The absolute errors for all 64 axially deformed nuclei in the mass regions 
A ^ 150 — 180 and A ^ 230 — 250 are smaller than 1 MeV. With stronger binding in 
symmetric nuclear matter (i.e. by increasing the absolute value of a^), the corresponding 
deviations of calculated binding energies become larger, and they also acquire a definite 
isotopic dependence (cf. Figs. IHl fTTIl . Reducing the absolute value of ay reverses the isotopic 
trend of the errors (cf. Fig. [6]). The isospin and mass dependence of binding energies shown 
in Figs. [6] - [H] is one of the main results of the present analysis, and it illustrates how 
sensitive are the calculated masses to the choice of the nuclear matter binding energy at sat- 
uration. It also clearly shows why it is not possible to accurately determine the parameters 
of a nuclear energy density functional already at nuclear matter level, i.e. in an ab initio 
approach starting from empirical NN and NNN interactions, without additional adjustment 
to low-energy data on finite medium-heavy and heavy nuclei. 

The results of Figs. [6]-[II]can be compared to those obtained with one of the most success- 
ful finite-range meson-exchange effective interactions: DD-ME2 (cf. Fig. [T2|) . The volume 
energy coefficient of DD-ME2 is ay = —16.14 MeV, and the parameters were adjusted to 
binding energies, charge radii and neutron radii of 12 spherical nuclei |30]- One notices both 
the pronounced isotopic and mass dependence of the deviations of binding energies calcu- 
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lated with DD-ME2. Although the span of the DD-ME2 absolute errors is somewhat smaller 
than the one of the corresponding point-coupling effective interaction with = —16.14 MeV 
(cf. Fig. [TTl) . the meson-exchange interaction obviously underbinds most of the 64 axially 
deformed nuclei, especially in the mass region A ^ 150 — 180. This is because DD-ME2 
was adjusted to binding energies of spherical nuclei, and therefore implicitly includes closed- 
shell effects beyond the mean field. Virtually all self-consistent relativistic models based on 
the static mean-field approximation are characterized by relatively small effective nucleon 
masses. The reason is that in the RMF framework the nonrelativistic-type effective mass 
^jvr 42] is not independent from the Dirac mass = m — S5. The latter determines 
not only the nucleon spin-orbit potential, but also the binding of symmetric nuclear matter 
and, therefore, constraints the nonrelativistic-type effective mass to a rather narrow interval 
around m*j^f^ ^ 0.65 m. A small effective mass translates into low density of single-nucleon 
states around the Fermi surface. This is especially pronounced in magic nuclei where stan- 
dard RMF models predict much too large energy gaps between occupied and unoccupied 
major shells. When these interactions are nevertheless forced to reproduce experimental 
binding energies of magic nuclei, i.e. when their parameters are adjusted to masses of nuclei 
like ^'^^Sn and ^''^Pb, the surrounding open-shell nuclei are predicted under-bound, giving 
rise to characteristic "arches" of the deviations between theoretical and experimental bind- 



mg energies 
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30|. "Arches" between shell closures, i.e. the over-binding of closed-shell 



nuclei relative to surrounding open-shell nuclei, characterize also most nonrelativistic self- 
consistent mean- field models, e.g. Skyrme-type effective interactions [15]. 

It has become customary to adjust nonrelativistic and relativistic energy density func- 
tionals to ground-state data of magic, closed-shell nuclei. However, the resulting effective 
interactions are not often used to calculate low-lying spectra of spherical nuclei. These 
functionals are more successful in the description of the evolution of deformation, shape 
coexistence phenomena, rotational bands, etc. in deformed nuclei. This is the rationale be- 
hind the present adjustment of the relativistic energy density functional directly to masses 
of axially deformed medium-heavy and heavy nuclei. This procedure, of course, does not 
solve the problem of "arches". They are still present, but now magic, closed-shell nuclei are 
over-bound with respect to experimental binding energies. 

In this section it has been shown that, among the effective density-dependent point- 
coupling interactions considered in the present analysis, the one with volume energy = 
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—16.06 MeV, surface energy = 17.498 MeV, and symmetry energy {S2) = 27.8 MeV 
(04 = 33 MeV), yields best results for the binding energies of axially deformed nuclei in the 
mass regions A ^ 150 — 180 and A ^ 230 — 250. We will denote this interaction DD-PCl 
(density-dependent point-coupling). In addition to the parameters of the isoscalar terms 
(set C in Table [T]), and the strength of the derivative term 6s = —0.815 fm^, DD-PCl 
is completely specified by the two parameters of the isovector channel: bxv = 1-836 fm^ 
and (Itv = 0.64. The total number of parameters is 10. In the next section a number of 
calculations will be performed that illustrate not only the predictive power of DD-PCl, but 
also problems in the calculation of masses of spherical nuclei. 

Finally, in Fig. [12] we compare the density dependence of the DD-PCl isoscalar vector 
and scalar nucleon self-energies in symmetric nuclear matter, with the starting approxima- 



;ion: the Hartree-Fock (HF) self-energies 33| calculated from the Idaho N^LO NN-potential 



3J] . As already emphasized above, the analysis of Ref . 33| has shown that at the relativis- 



tic Hartree-Fock level microscopic NN-potentials do not yield nuclear matter saturation. To 
achieve saturation of homogeneous symmetric matter, and to reproduce binding energies of 
finite nuclei, the self-consistent DD-PCl mean-fields must include effects of short-range cor- 
relations. This necessitates an increase in magnitude of the HF scalar self-energy for nucleon 
densities below 2psat- At saturation density, in particular, this increase is 70 MeV and we 
also note the pronounced exponential dependence on density of the DD-PCl self-energies, as 
compared to the almost linear density dependence of the HF Idaho self-energies. The mag- 
nitude of the scalar self-energy determines also the effective Dirac mass and, therefore, the 
strength of the effective nucleon spin-orbit potential in finite nuclei. At low densities below 
Psat the HF Idaho vector self-energy is much less modified by the requirement of saturation 
and self-consistency. An interesting result is that, at saturation density, the HF Idaho and 
DD-PCl vector self-energies differ by less than 3 MeV. At much higher nucleon densities the 
behavior of the DD-PCl self-energies has been determined by fixing the EoS to the point 
p = 0.56 fm~^ on the microscopic EoS of Akmal, Pandharipande and Ravenhall 36|, and 
therefore it can no longer be compared with the HF Idaho self-energies. 
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IV. ILLUSTRATIVE CALCULATIONS 



We have performed a series of test calculations of binding energies, charge isotope shifts, 
deformations, and isoscalar and isovector giant resonances. Ground-state properties are 
calculated using the axially deformed RMF model. Pairing correlations are treated in the 
BCS approximation with constant pairing gaps determined from the 5-point formula: 

A(5)(Aro) = [E{No + 2) - 4E{No + 1) + 6E{No) - 4E{No - 1) + E{No - 2)] , (18) 

8 

where E{Nq) denotes the experimental binding energy of a nucleus with A^o neutrons {Zq 
protons), and U^o = +1 (—1) for Nq even (odd). 



The relativistic quasiparticle random-phase approximation (RQRPA) [4J, |45|] is used 
to calculate excitation energies of giant resonances in spherical nuclei. Results calculated 
with DD-PCl are compared to available data, and with predictions of the meson-exchange 
interaction DD-ME2. 

In Fig. [H] the RMF-I-BCS predictions for charge radii of the Nd, Sm, Gd, Dy, Er and Yb 
isotopic chains are compared with data from Ref. 46|]. The charge density is obtained by 
folding the theoretical point-proton density distribution with the Gaussian proton-charge 
distribution. For the latter an rms radius of 0.8 fm is used, and the resulting ground-state 
charge radius reads 

rc = ^Jr^ + 0.64 fm , (19) 

where Vp is the radius of the point-proton density distribution. Even though the two RMF 
interactions, meson-exchange DD-ME2 and point-coupling DD-PCl, have been adjusted 
using different procedures and to different data sets, they predict virtually identical charge 
radii for all six isotope chains. The theoretical values are in excellent agreement with data 
for Nd, Sm, Gd, and Dy nuclei. For the heavier isotopes of Er and Yb, the calculated 
radii are only slightly above the experimental values. Note that the parameters of DD- 
ME2 were tuned both to binding energies and charge radii of spherical nuclei, whereas only 
experimental masses of deformed nuclei have been used to adjust the interaction DD-PCl. 
Of the nuclides that belong to the six isotopic chains shown in Fig. [TU only those with 
> 92 are included in the data set of 64 deformed nuclei used to determine the density 
functional DD-PCl. 

The ground-state quadrupole deformation parameters P2 are calculated according to the 
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prescription of Ref. [47|. The theoretical predictions for the quadrupole deformation pa- 
rameters of Nd, Sm, Gd, Dy, Er and Yb isotopes are displayed in Fig. [151 iii comparison 
with the empirical values from Ref. 48|]. Also in this case DD-ME2 and DD-PCl predict 
identical ground-state shapes for all six isotopic chains and reproduce not only the global 
trend of the data, but also the saturation of quadrupole deformation for heavier isotopes. 
The only discrepancy is found around the = 82 closed shell in Nd and Sm isotopes, where 
both interactions predict spherical ground states, whereas data indicate that these nuclei are 
slightly prolate deformed, probably with soft potential energy surfaces. Shape coexistence 
structures (spherical-deformed, or prolate-oblate shapes), and soft potential surfaces in gen- 
eral, cannot quantitatively be described on a mean-field level. The description of coexisting 
shapes must include long-range correlations, and necessitates angular momentum projec- 
tion and configuration mixing [is], [l4|, not considered in this work. For heavier Nd and 
Sm isotopes, however, the predictions of both DD-ME2 and DD-PCl are again in excellent 
agreement with empirical prolate deformations. 

Even though DD-PCl is not constructed with the idea of being used as a mass formula, 
nevertheless this functional must also be tested in the calculation of binding energies. We 
consider the cases of deformed and spherical nuclei separately. As a first test Fig. [16] shows 
the absolute deviations of the DD-PCl binding energies from experimental values of de- 
formed nuclei in the mass regions A ^ 120 — 130, A ^ 150 — 180 and A ^ 230 — 250, as 
functions of the asymmetry coefficient and mass number. Cross symbols denote the 64 nuclei 
used to adjust the parameters of DD-PCl, and correspond to the deviations already shown 
in Fig. [71 Additional deformed nuclei that have not been used in the fit, are represented 
by diamond symbols. We include about 20 nuclei in the mass region A ^ 120 — 130, and 
12 more around mass A ^ 150 — 160. In this way the predictions for binding energies are 
extrapolated to a lower mass region not included in the fit and, even more importantly, to 
lower values of the asymmetry parameter. The overall agreement with data is very good, 
and the absolute deviations from experiment are contained in the interval ±1 MeV. 

In the case of spherical closed-shell nuclei, the variance between calculated masses and the 
corresponding experimental values is somewhat larger. This is illustrated in Fig. [T71 where, 
for the Pb and Sn isotopic chains, we plot the absolute deviations of the calculated binding 
energies from data as functions of the mass number. The binding energies calculated using 
the RMF-I-BCS model with the functional DD-PCl are also compared to those obtained 
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with the meson-exchange interaction DD-ME2. The latter, hke most modern self-consistent 
mean-field nonrelativistic and relativistic interactions, was adjusted to reproduce the binding 
energies of doubly closed-shell nuclei, including ^'^^Sn and ^°'^Pb. In addition, of the isotopes 
shown in Fig. [TTl the set of nuclei used to fine-tune the parameters of DD-ME2 included 
also ^^^Sn, ^^^Sn, ^°^Pb and ^^^Pb. The resulting binding energies of the two isotopic chains 
are, of course, in better agreement with data than those predicted by the density functional 
DD-PCl. Because it has been tailored to masses of deformed nuclei, DD-PCl necessarily 
overbinds spherical closed-shell nuclei. The situation is actually not as bad as the comparison 
with the experimental mass of ^^^Sn might indicate. This isotope, with a deviation of 5.21 
MeV, is in fact the worst case of those nuclei that we have calculated so far. For instance, ^^O 
is calculated overbound by 0.72 MeV, ^^Ca by 1.51 MeV, ^o^Pb by 3.51 MeV, etc. Although 
the interaction DD-ME2 predicts masses of spherical nuclei closer to data, it underbinds 
most deformed nuclei (cf. Fig. [T2!). 

The origin of the additional binding in closed-shell nuclei can be found in the predicted 
shell structure. For the two particular examples of ^osp]^ g^,^^ ^^^Sn, this effect is illustrated 
in Fig. [18], which shows the comparison between experimental and DD-PCl single-nucleon 
spectra of protons (upper panel) and neutrons (lower panel). The experimental spectra 
correspond to the single-nucleon separation energies of Ref. |49|. Note that single-nucleon 
orbitals are solutions of the relativistic Kohn-Sham equations and the corresponding eigen- 
values, introduced just as Lagrange multipliers, do not have a directly observable physical 
interpretation, i.e. they cannot exactly be identified with nucleon separation energies 50|. 
For the last few occupied orbitals close to the Fermi surface, however, the Kohn-Sham 
eigenvalues approximately correspond to physical single-nucleon energies. As with most 
self-consistent mean-field models jl|, the calculation reproduces the overall structure and 
ordering of single-nucleon levels, but not the level density around major shell gaps. Be- 
cause of the low effective nucleon mass, the magnitude of the spherical shell gaps between 
occupied and unoccupied states is overestimated. As Fig. [18] clearly shows, the theoretical 
occupied levels are on the average considerably deeper than the corresponding empirical 
single-nucleon states, and this effect gives rise to the overbinding that characterizes masses 
of spherical nuclei calculated with nonrelativistic and relativistic mean-field models. We 
note that a similar analysis for the nonrelativistic Skyrme interaction Sly4 was carried out 
in Ref. [15 1. 
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This effect is particularly pronounced for the = 82 neutron gap in ^^^Sn. Not only 
is the theoretical gap much larger than the empirical one, but the functional DD-PCl, like 
many other mean-field interactions, does not reproduce the empirical sequence of neutron 
levels below = 82. The last occupied level should be 2d3/2, whereas DD-PCl places this 
orbital below lhii/2 and 3si/2- This increases the magnitude of the gap, and the net result 
is the additional binding of heavy Sn isotopes shown in Fig. [T71 If an interaction with a low 
effective nucleon mass, like DD-PCl or DD-ME2, is nevertheless tuned to the masses of ^^^Sn 
and ^''^Pb, most deformed nuclei in between will be underbound, as shown for DD-ME2 in 
Fig. [El 

Even though it will not be considered in the present work, we would like to point out 
that to enhance the nonrelativistic effective (Landau) nucleon mass, the functional must 
go beyond the static mean-field approximation, and include momentum- dependent (energy- 
dependent in stationary systems) nucleon self-energies. In the meson-exchange RMF frame- 
work, in particular, this has been achieved by including in the effective Lagrangian a particu- 



lar form of coupling between meson fields and the derivatives of the nucleon fields 5l|, l52l, l53| 



This leads to a linear momentum dependence of the scalar and vector self-energies in the 
Dirac equation for the in-medium nucleon. Although this extension of the standard mean- 
field framework is phenomenological, it is nevertheless based on Dirac-Brueckner calculations 
of in-medium nucleon self-energies, and consistent with the relativistic optical potential in 
nuclear matter, extracted from elastic proton-nucleus scattering data. In this way it was 
possible to increase the effective (Landau) mass to ^jvi? ~ 0.8 m. An additional enhance- 
ment of the effective nucleon mass in finite nuclei is caused by the coupling of single-nucleon 
levels to low-energy collective vibrational states 5J|, an effect which goes entirely beyond the 
mean-field approximation. In the RMF framework the coupling of single-nucleon states to 
low-energy phonons, and the resulting increase of the effective mass, were recently explored 
in Ref. [55 1. 

The somewhat more pronounced deviations between the theoretical DD-PCl and experi- 
mental masses of the Sn and Pb isotopic chains, do not affect the accuracy of the calculated 
radii of these nuclei. In Fig. [19] the RMF+BCS model predictions for the charge radii of Pb 
and Sn isotopes, calculated with the effective interactions DD-PCl and DD-ME2, are shown 
in comparison with empirical values [46^. The two interactions predict virtually identical 
values for the charge radii, in excellent agreement with data for the Sn nuclei, and only 



21 



slightly above the empirical charge radii of the Pb isotopes. We note again that the charge 
radii of ^^^Sn, ^^^Sn, ^^^Sn, ^°^Pb, ^''^Pb, and ^^^Pb were used to tune the parameters of 
DD-ME2, whereas DD-PCl has only been adjusted to masses. A similar result is also ob- 
tained for the thickness of the neutron skin in Pb and Sn nuclei. In Fig. [2D] the DD-ME2 and 
DD-PCl results for the differences between the neutron and proton rms radii are compared 



with available data 
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581]. Although the experimental uncertainties are large, both 



interactions nicely reproduce the isotopic trend of neutron radii in Sn nuclei, and predict 
values in very good agreement with the empirical values of neutron skin thickness. These 
data, however, were specifically used to adjust the isovector channel of the DD-ME2 inter- 
action, whereas in the case of DD-PCl this level of agreement is achieved with the choice of 
the asymmetry energy at saturation 04 = 33 MeV. 

Another very important test of self-consistent mean-field models are excitation energies 
of collective modes and, in particular, giant multipole resonances. Using the relativistic 
(Q)RPA 4J, |45l] with the DD-PCl functional, we have therefore carried out few representa- 
tive calculations of giant resonances in spherical nuclei. The RQRPA is fully self-consistent, 
i.e. both in the particle-hole and particle-particle channels the same interactions are used 
in the equations that determine the ground state of a nucleus, and as residual interac- 
tions in the matrix equations of RQRPA. The RQRPA configuration space includes also 
the Dirac sea of negative energy states. In adjusting the parameters of DD-PCl we took 
into account the results of Refs. 35|], where it has been shown that, to reproduce the ex- 
perimental excitation energies of isoscalar giant monopole resonances (ISGMR) in spherical 
nuclei, relativistic point-coupling interactions require a nuclear matter compression modu- 
lus of Knm ~ 230 MeV, somewhat lower than the values typically used for meson-exchange 
relativistic interactions J39I, and within the range of values used by modern nonrelativis- 
tic Skyrme interactions. In Ref. 39| it was also shown that data on the isovector giant 
dipole resonance (IVGDR) constrain the range of the nuclear matter symmetry energy at 
saturation to 31 MeV < 04 < 35 MeV, and in this work we have used 04 = 33 MeV for 
all point-coupling interactions. For ^°'^Pb the RRPA results for the isoscalar monopole and 
isovector dipole response are shown in Fig. [21] For the multipole operator Qx^ the response 
function R{E) is defined 
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where |0) denotes the ground state of an even-even nucleus, |Aj) is the corresponding i-th 
excited state of multipolarity A, F is the width of the Lorentzian distribution, and 

B{X^O) = ^^m\Qx\\X)\'. (21) 
The k-th moment of the strength function is defined by: 

m,{Q,) = J2 El \{MQm^ (22) 

i 

The ratio mi/rriQ is the quantity often compared with the experimental excitation energy 
of the corresponding resonance although, of course, this is strictly correct only if there 
are no pronounced multiple peaks within the energy interval over which the sumation in 
Eq. fl22l) is performed. In the examples considered here the continuous strength distributions 
are obtained by folding the discrete spectrum of R(Q)RPA states with the Lorentzian (cf. 
Eq. ([20])) of constant width T = 1 MeV. 

The calculated excitation energies of the giant resonances in ^osp]^ ^.g^j^ compared with 
very accurate data. For the ISGMR the calculated rru/rriQ = 14.17 MeV is rather close to 
the experimental value mi/niQ = 13.96 ± 0.2 MeV 59|]. The relativistic RPA peak energy 
of the IVGDR at 13.6 MeV is also in very good agreement with the experimental excitation 



energy E* = 13.3±0.1 MeV [60(. A similar level of agreement, both for ISGMR and IVGDR, 
is also obtained with the DD-ME2 interaction, however we note that these data were taken 
into account in adjusting the parameters of DD-ME2 jsol]. 

In Fig. [22] the RQRPA results for the Sn isotopes are compared with data on IVGDR 
excitation energies |61|]- In contrast to the case of ^^'^Pb, the strength distributions in the 
region of giant resonances exhibit considerable fragmentation. Note, however, that the 
RQRPA calculation with the DD-PCl interaction reproduces the isotopic dependence of 
the IVGDR and the experimental excitation energies. The theoretical peak energies 15.56 
MeV (^^^Sn), 15.35 MeV ("^Sn), 15.26 MeV {^^°Sn), 15.13 MeV {^^^Sn), are in excellent 
agreement with the experimental values: 15.68 MeV (^^^Sn), 15.59 MeV (^^^Sn), 15.36 MeV 
(^20Sn), 15.19 MeV (^^^Sn), respectively. 

Finally, in Fig. [23] we display the RQRPA isoscalar monopole strength functions for the 
chain of even-even Sn isotopes: ^^^~^^^Sn. The evolution of ISGMR in Sn isotopes can be 



compared with very recent data from Ref. 62]. In general, the theoretical excitation energies 



Eq = mi/ mo, when compared with the corresponding experimental values evaluated in 
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the energy interval between 10.5 MeV and 20.5 MeV, are systematically somewhat higher: 
between 0.8 MeV and 1 MeV. This result might indicate that the value of the nuclear matter 
compression modulus for this functional could actually be chosen Knm < 230 MeV. However, 

Sn isotopes are much lighter than ^°^Pb, and one can expect that the calculated ISGMR will 
be more affected by the surface incompressibility, a quantity that we have not attempted to 
determine in this work. 

V. SUMMARY AND CONCLUDING REMARKS 

The principal goal of modern nuclear structure theory is to build a consistent microscopic 
framework that will describe ground-state properties, nuclear excitations and reactions at 
a level of accuracy comparable with experimental results, and provide reliable predictions 
for systems very far from stability, including data for astrophysical applications that are 
not accessible in experiments. At present the only viable approach to a comprehensive 
description of arbitrarily heavy nuclear systems, including vast regions of short-lived nuclei 
with extreme isospin values and extended nucleonic matter, is the one based on nuclear 
energy density functionals. 

In this work we have explored a particular class of relativistic nuclear energy density 
functionals in which only nucleon degrees of freedom are explicitly used in the construction 
of effective interaction terms. Short-distance correlations, as well as intermediate and long- 
range dynamics, are effectively taken into account in the nucleon density dependence of the 
strength functionals of second-order contact interactions in an effective Lagrangian. Start- 
ing from microscopic nucleon self-energies in nuclear matter, a phenomenological ansatz for 
the density-dependent coupling functionals has been formulated and the corresponding pa- 
rameters adjusted in self-consistent mean-field calculations of masses of 64 axially deformed 
nuclei in the mass regions A fa 150 — 180 and A ^ 230 — 250. 

The relationship between global properties of the nuclear matter equation of state (EoS) 
and the corresponding predictions for nuclear masses has been analyzed in detail. Ground 
states of deformed nuclei have been calculated in the self-consistent mean field approximation 
by employing sets of effective interactions with different values of the volume a^, surface a^, 
and symmetry energy 0-4 in nuclear matter, whereas empirical constraints have been used 
for the nuclear matter saturation density, compression modulus, and Dirac effective mass. 
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The calculated masses are not particularly sensitive to the saturation density, and previous 
relativistic RPA calculations of excitation energies of isoscalar giant monopole and isovector 
giant dipole resonances in finite nuclei, as well as results for the neutron skin thickness have 
been used to determine the nuclear matter compression modulus and symmetry energy, 
respectively 

One of the important results of the careful analysis of deviations between calculated and 
experimental masses (mass residuals), is the pronounced isospin and mass dependence of 
the residuals on the nuclear matter volume energy at saturation. To reduce the absolute 
mass residuals to less than 1 MeV, and to contain their mass and isotopic dependence, strict 
constraints on the value of must be met. The narrow window of allowed values of the 
volume energy cannot be determined microscopically already at the nuclear matter level, 
but rather results from a fine-tunning of the parameters of the energy density functional to 
experimental masses. Calculated binding energies and charge radii are also sensitive to the 
choice of the surface coefficient that determines the surface energy and surface thickness 
of semi-infinite nuclear matter. In the functional considered in the present work, these quan- 
tities are controlled by the strength of the derivative isoscalar-scalar coupling interaction. 
The range of allowed values of the strength parameter determined by data is in very good 
agreement with estimates obtained from microscopic calculations of inhomogeneous nuclear 
matter. 

The optimal energy density functional (DD-PCl) determined in a multistep parameter 
fit to the masses of 64 axially deformed nuclei, has been further tested in calculations of 
properties of spherical and deformed medium-heavy and heavy nuclei, including binding en- 
ergies, charge radii, deformation parameters, neutron skin thickness, and excitation energies 
of giant multipole resonances. Results have been compared with available data, and with 
predictions of one of the most successful finite-range meson-exchange relativistic effective in- 
teractions: DD-ME2. In general, a very good agreement with data has been obtained except, 
perhaps, for the effect of overbinding of spherical closed-shell nuclei. DD-PCl, like virtu- 
ally all relativistic mean-field models, is characterized by a relatively low effective nucleon 
mass and, when adjusted to masses of deformed nuclei, it overbinds spherical closed-shell 
systems. The well known problem of "arches" of mass residuals between shell closures could 
be addressed by a functional that goes beyond the static mean-field approximation, but 
this approach has not been considered in the present model. Very good results have been 
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obtained for the excitation energies of giant monopole and dipole resonances in spherical 
nuclei, calculated with the relativistic quasiparticle random-phase approximation based on 
the DD-PCl functional. The agreement with data confirms the choice of the nuclear matter 
compressibility and symmetry energy for DD-PCl. 

The total number of parameters in the functional DD-PCl is 10, like in most non- 
relativistic Skyrme-type density functionals. Note, however, that the effective Lagrangian 
of DD-PCl contains only four interaction terms except, of course, the Coulomb term (cf. 
Eq. ([2])). The 10 parameters determine the density dependence of the strength functionals 
and reflect the complex nuclear many-body dynamics. We also emphasize that, because 
the high-density behavior of the corresponding nuclear matter EoS has been adjusted to a 
microscopic equation of state extensively used in studies of high-density nucleon matter and 
neutron stars, DD-PCl should also be tested in astrophysical applications. 

This work is part of a program to develop an universal relativistic energy density func- 
tional to be used in studies of of the evolution of shell structure, deformation, shape co- 
existence phenomena, rotational bands, etc. in transitional and deformed nuclei. In the 
first step the parameters of the density functional have been adjusted to reproduce bind- 
ing energies of a large set of axially deformed nuclei. In the continuation of this program 
we plan to build a functional, based on DD-PCl, to be used in the new relativistic model 
that employs the generator coordinate method (GCM) to perform configuration mixing of 
angular-momentum and particle-number projected mean-field wave functions. However, if 
rotational energy corrections and quadrupole fluctuations are treated explicitly in the GCM 
framework, they should not at the same time implicitly be included in the functional, i.e. 
through parameters adjusted to data that already include correlations. Therefore, starting 
from DD-PCl, the parameters of this new functional can be adjusted to pseudodata, ob- 
tained by subtracting correlation effects from experimental masses and radii. The resulting 
energy density functional will be tested in relativistic GCM model studies of shell evolution, 
deformations, shape coexistence and shape phase transitions. 
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FIG. 1: (Color online) The equations of state of symmetric nuclear matter (binding energy as a 
function of nucleon density) for the eight point-coupling effective interactions of Table HI in compar- 



ison with the EoS of the meson-exchange effective interaction DD-ME2 



EoS of Ref. 



30(], and the microscopic 



36l |. The two points from the microscopic EoS on which the point-coupling effective 



interactions (sets A-H) were adjusted, are denoted by larger filled circle symbols. 
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FIG. 2: (Color online) Surface energy of semi-infinite nuclear matter as a function of the surface 
thickness, for the eight point-coupling effective interactions of Table [H The corresponding values 
of the strength Ss of the derivative coupling term in the point-coupling Lagrangian Eq. ([2]), are 
displayed on the upper horizontal axis. The filled square symbol denotes the surface energy and 



surface thickness predicted by the meson-exchange effective interaction DD-ME2 
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FIG. 3: (Color online) Surface energies of semi-infinite nuclear matter that minimize the deviation 
of the calculated binding energies from data, for the set of nuclei listed in Table [Til plotted as 
functions of the volume energy at saturation, for six values of the symmetry energy (5*2 ). 
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FIG. 4: (Color online) x^-deviations Eq. ()17p of the theoretical binding energies from data for 
the set of deformed nuclei listed in Table |lTl and for each combination of the surface, volume and 
symmetry energy shown in Fig. [3j 



33 



160 


1 1 


1 1 1 1 


140 
120 
100 


<S2>=27.6 MeV 
<S2>=27.8 MeV 
<S2>=28.0 MeV 
<S2>=28.2 MeV 

<S2>=28.6 MeV 


■ / 

• / 

T / 

♦ / 


80 






60 






40 






20 






n 


1 1 


1 1 1 1 



16.04 16.06 16.08 16.10 16.12 16.14 
|aj (MeV) 



FIG. 5: (Color online) The minimum x'^-deviation of the theoretical binding energies from data, as a 
function of the volume energy coefficient. Each point represents the minimum of the corresponding 
curve plotted in Fig. [H 
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FIG. 6: (Color online) Absolute deviations of the calculated binding energies from the experimental 
values of the 64 axially deformed nuclei listed in Table HH as functions of the asymmetry coefficient 
(upper panel), and mass number (lower panel). Lines connect nuclei that belong to the isotopic 
chains shown in the legend. The theoretical binding energies are calculated using the point-coupling 
effective interaction characterized by the volume energy = —16.04 MeV. 
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FIG. 7: (Color online) Same as in Fig. [H for the point-coupling effective interaction with volume 
energy = —16.06 MeV. 
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FIG. 8: (Color online) Same as in Fig. [H for the point-coupling effective interaction with volume 
energy = —16.08 MeV. 
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FIG. 9: (Color online) Same as in Fig. [H for the point-coupling effective interaction with volume 
energy = —16.10 MeV. 
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FIG. 10: (Color online) Same as in Fig. [6l for the point-coupling effective interaction with volume 
energy = —16.12 MeV. 
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FIG. 11: (Color online) Same as in Fig. [6l for the point-coupling effective interaction with volume 
energy = —16.14 MeV. 
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FIG. 12: (Color online) Absolute deviations of the calculated binding energies from the experi- 
mental values of the 64 axially deformed nuclei listed in Table [III as functions of the asymmetry 
coefficient (upper panel), and mass number (lower panel). Lines connect nuclei that belong to 
the isotopic chains shown in the legend. The theoretical binding energies are calculated using the 
meson-exchange effective interaction DD-ME2 3C|. 
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FIG. 13: (Color online) Vector (upper panel) and scalar (lower panel) nucleon self-energies in 
symmetric nuclear matter as functions of the nucleon density. The self-energies that correspond to 
the phenomenological density functional DD-PCl are compared with the Hartree-Fock self-energies 



33| calculated from the Idaho N'^LO NN-potential 
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FIG. 14: (Color online) Charge radii of Nd, Sm, Gd, Dy, Er and Yb isotopic chains. The results of 
the RMF+BCS calculation with the DD-PCl and DD-ME2 interactions are compared with data 
46|. 
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FIG. 15: (Color online) DD-PCl and DD-ME2 predictions for the ground-state quadrupole de- 
brmations /?2 of the Nd, Sm, Gd, Dy, Er and Yb isotopes, in comparison with empirical values 
4£|. 
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FIG. 16: (Color online) Absolute deviations of the DD-PCl binding energies from experimental 
values of deformed nuclei in the mass regions A ^ 120 — 130, A 150 — 180 and A 230 — 250, 
as functions of the asymmetry coefficient (upper panel), and mass number (lower panel). 
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FIG. 17: (Color online) Absolute deviations of the calculated binding energies from experimental 
values for the Pb (upper panel) and Sn (lower panel) isotopic chains, as functions of the mass 
number. The theoretical binding energies are calculated using the RMF+BCS model with the 
point-coupling effective interaction DD-PCl, and the finite-range meson exchange interaction DD- 
ME2. 
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FIG. 18: (Color online) Comparison between experimental (left) and DD-PCl (right) single-nucleon 
spectra of protons (upper panel) and neutrons (lower panel), for ^ospb and ^^^Sn. The experimental 
spectra are from Ref. [49]. 
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FIG. 19: (Color online) RMF+BCS model predictions for the charge radii of Pb (upper panel) and 
Sn (lower panel) isotopes, calculated with the DD-PCl and DD-ME2 effective interactions, and 
compared with data j46| . 
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FIG. 20: (Color online) RMF+BCS model predictions for the differences between the neutron and 



proton rms radii of Pb (upper panel) and Sn (lower panel) isotopes, calcu 



ated with the DD-PCl 



and DD-ME2 effective interactions, in comparison with available data 
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FIG. 21: (Color online) The isoscalar monopole (left panel) and isovector dipole (right panel) 
strength distribution in ^''^Pb calculated with the relativistic RPA using the effective interaction 
DD-PCl. The experimental excitation energy are denoted by arrows: 13.96 it 0.2 [59,] for the giant 
monopole resonance, and 13.3 ±0.1 for the giant dipole resonance, respectively. 
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FIG. 22: (Color online) The RQRPA isovector dipole strength functions in Ii6,ii8,i20,i24g^^ calcu- 
lated with the DDPCl effective interaction. The experimental IVGDR excitation energies joil are 
denoted by arrows. 
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FIG. 23: (Color online) The RQRPA isoscalar monopole strength distributions in even-even 
ii2-i24gj^ nuclei, calculated with the DDPCl effective interaction. Arrows denote the positions 
of experimental ISGMR excitation energies 62 1. 
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TABLE I: Isoscalar parameters Eq. (jl3p of the point-coupling effective interactions with volume 
energy coefficients: = -16.02 MeV (set A), ay = -16.04 MeV (set B), Uy = -16.06 MeV (set 
C),ay = -16.08 MeV (set D) , a^, = -16.10 MeV (set E), a„ = -16.12 MeV (set F), a„ = -16.14 
MeV (set G) and ay = -16.16 MeV (set H). 



parameter 


set A 


set JD 


oof 
set kj 


set JJ 


set Hj 


set r 


set (j 


set rl 


as (fm2) 


-10.0220 


-10.0332 


-10.0462 


-10.0855 


-10.0951 


-10.1051 


-10.1137 


-10.1220 


bs (fm2) 


-9.1781 


-9.1666 


-9.1504 


-9.0623 


-9.0539 


-9.0436 


-9.0384 


-9.0307 


cs (fm^) 


-6.2799 


-6.3541 


-6.4273 


-6.4878 


-6.5611 


-6.6336 


-6.7065 


-6.7786 


ds 


1.3585 


1.3654 


1.3724 


1.3806 


1.3872 


1.3938 


1.4001 


1.4065 


ay (fm^) 


5.9020 


5.9108 


5.9195 


5.9262 


5.9348 


5.9431 


5.9513 


5.9594 


by (fm2) 


8.8711 


8.8687 


8.8637 


8.8156 


8.8150 


8.8134 


8.8148 


8.8147 


Cv (fm^) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


dv 


0.6548 


0.65676 


0.6584 


0.6547 


0.6568 


0.6587 


0.6610 


0.6630 



TABLE IL The binding energies of the isotopic chains 62 < Z < 72 and 90 < Z < 98 have been 
used to adjust the parameters of relativistic point-coupling effective interactions. Nmin and Nmax 
denote the corresponding ranges of neutron number in even-even nuclides. 

Z 62 64 66 68 70 72 90 92 94 96 98 
N.min 92 92 92 92 92 72 140 138 138 142 144 
Nmax 96 98 102 104 108 110 144 148 150 152 152 
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